Compositional uniformity, domain patterning and the mechanism underlying 

nano-chessboard arrays 

Santiago Gonzalez 1 , J.M. Perez-Mato 1 ^] Luis Elcoro 1 , Alberto Garcia 2 , Ray L. Withers 3 , and L. Bourgeois 4 
1 Departamento de Fisica de la Materia Condensada, Facultad de Ciencia y Tecnologia, 
Universidad del Pais Vasco (UPV-EHU), Apdo. 644, 48080 Bilbao, Spain. 
2 Institut de Ciencia de Materials de Barcelona, 
(ICMAB-CSIC), Campus de la UAB, E-08193 Bellaterra, Spain. 
3 Research School of Chemistry, The Australian National University, Canberra, A.C.T, Australia and 
4 Monash Centre for Electron Microscopy (MCEM), 
Monash University, Clayton, Victoria, Australia. 
(Dated: July 12, 2012) 

We propose that systems exhibiting compositional patterning at the nanoscale, so far assumed to 
be due to some kind of ordered phase segregation, can be understood instead in terms of coherent, 
single phase ordering of minority motifs, caused by some constrained drive for uniformity. The 
essential features of this type of arrangements can be reproduced using a superspace construction 
typical of uniformity-driven orderings, which only requires the knowledge of the modulation vectors 
observed in the diffraction patterns. The idea is discussed in terms of a simple two dimensional 
lattice-gas model that simulates a binary system in which the dilution of the minority component is 
favored. This simple model already exhibits a hierarchy of arrangements similar to the experimen- 
tally observed nano-chessboard and nano-diamond patterns, which are described as occupational 
modulated structures with two independent modulation wave vectors and simple step-like occupa- 
tion modulation functions. 

PACS numbers: 61.50.Nw, 61.44.Fw, 62.23.Pq,68.37.Lp 



I. INTRODUCTION 

Domain self-patterning on the nanoscale, particularly 
if the nanoscale chemical ordering can be engineered or 
tuned, is the goal of much research in the nanoscience 
and nanotechnology area. Potential applications include 
the use of the surfaces of such materials as templates 
for the assembly of molecular monolayers, for the re- 
liable synthesis of functional nano-structured materials 
or for the structured adsorption of gas species. The 
recent discovery that tuneable nanoscale ordering oc- 
curs spontaneously in the wide range, non-stoichiomctric 
Lii/2-3xLn 1 / 2 + 2 ;Ti03, 0.02< x < 0.12, solid solution^ as 
well as in a variety of other A-sitc ordered perovskitc- 
relatcd phases — — is thus of intense current research 
interest. Guiton and Daviesi demonstrated the exis- 
tence of 2D, nanoscale compositional ordering in the 
Lii/ 2 -3 :c Ndi/2 +a; Ti03, 0.02 < x < 0.12, system via High 
Angle Annular Dark Field (HAADF) imaging and in- 
terpreted this result in terms of phase segregation aris- 
ing from some kind of "ordered" spinodal decomposition 
into Li!/2Ndi/2Ti03 nano-chessboard regions separated 
by narrow Nd2/3Ti03 boundary regions. As an exam- 
ple, Fig. [1] shows the kind of nano-patterns that can be 
observed in this compound. 

The 2D apparent phase separation into composition- 
ally distinct areas described above is strongly reminiscent 
of the ID atomic order that takes place in many compo- 
sitionally flexible ordered phases. The list is very exten- 
sive, but one can cite for instance many binary and/or 
ternary alloy systems, —— reduced rutile structures such 
as the Ti„02n-i system^i Re03 related compounds 



such as the MoOs-^ system^ the hexagonal perovskite 
related LaTii_ x 3 family^, the (1 - a;)Ta 2 5 + xW0 3 
system ) 10 ' 11 and others^ 2 - In most of these ID composi- 
tionally ordered arrangements, the basic driving mech- 
anism is the maximization of the separation of minor- 
ity structural motifs within a matrix of an underlying 
sub-structure. In the case of binary alloy systems, for 
example, there is usually a well-defined majority motif 
such as an ordered Cu 3 Au unit which is, as regularly as 
possible, interleaved with a minority anti-phase bound- 
ary (APB) type structural motif, which locally has a dif- 
ferent composition and leads to a continuously variable 
overall composition dependent upon the spatial distribu- 
tion of the minority motif. Likewise, in the case of the 
reduced rutile Ti„02n-i system, the majority motif is 
the rutile type Ti02 structure, while the minority mo- 
tif is a crystallographic shear plane (CSP), which again 
changes the overall composition if "regularly" and uni- 
formly arranged. These compositional orderings are one- 
dimensional and often yield patterns with composition 
varying in the nanoscale in the form of stripes (e.g. sec 

Fig. [5]). 

In the last decade it has been shown2iI2 that this type 
of one-dimensional uniformity-driven arrangement can be 
most efficiently described by means of step-like atomic 
occupation modulations with wave vectors and forms for 
the modulation functions that can be systematized ac- 
cording to simple rules using the superspace formalism. 
— _ — Very recently, the main ideas have been extended 
to the analysis of the distribution of minority motifs in 
two-dimensional systems^ showing that pseudouniform 
arrangements can originate as the result of a simple quest 
for maximal dilution of a minority motif within lattice- 
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FIG. 1: (a) A typical low resolution TEM image 

of the chessboard type nano-structure characteristic of a 
Lii/2-3a;Nd 1/ /2+ : ETi03, £=0.067, sample, (b) a high resolu- 
tion Bright Field STEM image of this nano-chessboard array. 
The lighter regions are Nd-poor and presumed to correspond 
to regions of stoichiometry Li 1 / 2 Ndi/2Ti03 while the darker 
boundary regions are Nd-rich and presumed to correspond to 
the stoichiometry Nd2/3Ti03, according to Ref. [J 

restricted sites. In wide ranges of composition, the dilu- 
tion drive was shown to stabilize ID stripe arrangements 
oblique to the underlying lattice, following predictable 
rules based on the alternation of sub-units of different 
compositions on the nanoscale, but still keeping an over- 
all identity as one single phase. The stripe arrangements 
were again shown to be easily describable as occupation- 
ally modulated structures, needing only a single modula- 
tion vector (as befits their intrinsic one-dimensionality) 
that can be determined unambiguously from the diffrac- 
tion diagram of the system. 

In this paper we show that these findings can be ex- 
tended to the case of intrinsically two-dimensional pat- 
terns, which need two modulation vectors for their su- 
perspace description. These 2D pseudouniform arrange- 
ments appear in certain ranges of composition, and in 
some cases exhibit nano-chessboard-like features simi- 
lar to those observed in the compounds mentioned at 
the beginning of this section. These patterns, resulting 
again from a basic drive to uniformity or maximal dilu- 
tion within a lattice, involve a coherent arrangement of 




FIG. 2: HRTEM image of a striped ordering in 

Bii_ a: Ca a; Fe03„ a; /2, ~ 0.20 < x <~ 0.49, a perovskite-related 
phase viewed along an [001] projection direction. Stripes of 
width 5a p are realized for this particular grain (for details see 
Ref.QJ). 

structural motifs, which is rather different from the idea 
of ordered spinodal decomposition or phase segregation 
that has been previously proposed. 

The paper is organized as follows. In Sec.|TI]we review 
the basic ideas and results on pseudouniform stripe order- 
ings in a simple two-dimensional A\_ X B X binary model 
presented recently (Ref. ITrh . In Sec. Mil we discuss the 
features and supcrspacc description of a different fam- 
ily of pseudouniform A\_ X B X orderings which are intrin- 
sically two-dimensional, including the determination of 
their two modulation vectors and their direct correspon- 
dence with the composition of the system. Finally, in 
Sec. IIVI we show how the above ideas can be useful to 
describe orderings in real systems. 



II. PSEUDOUNIFORM ID STRIPE 
COMPOSITIONAL ORDERING 

The basic patterns appearing due to a drive to max- 
imal dilution of minority motifs, abstracted from any 
other effects present in real systems, can be studied us- 
ing a simple model of a square lattice with binary com- 
position Ai_ x B x , in which the B "atoms" represent the 
minority motifs. In Ref. [13 the dilution drive was orig- 
inally mimicked by isotropic Yukawa-type repulsive in- 
teractions between the B particles, and a simple lattice- 
gas simulated-annealing method was used to explore and 
generate pseudouniform patterns by minimizing the re- 
pulsion energy. For most compositions, stripe orderings 
of the type shown on Fig. [3] were found to optimize the 
lattice-ordered dilution of the B atoms. The stripes are 
formed by the concatenation of tiles, parallelepipeds de- 
lineated by B atoms, with a fixed composition and ori- 
entation within a stripe. For example, in Figure [3ja), 
the B atoms (the larger black dots) are present in the 
proportion a;=5/12 and they achieve an ordered distri- 
bution as close as possible to uniformity by forming a 
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sequence of 73-rich and 73-poor stripes of tiles with com- 
position 1/2 and 1/3 respectively (the smaller, square 
tile contains one A and one 73 atom, while the larger 
tile contains one 73 and two A atoms). As the overall 
B fraction £=5/12 is closer to 1/2 than to 1/3, the 1/2 
stripe regions are dominant (by a ratio 3:2), and the 1/3 
stripes have minimal width and are isolated, as if they 
could be considered minority motifs themselves, subject 
to a dilution drive. The fact that the stripe sequence and 
makeup follow strict rules as a function of overall compo- 
sition was one of the main results of Rcf. The rules 
are based on the decomposition of the overall B frac- 
tion x in irreducible terms according to the Farey tree 
constructio n 17 ! 18 : for example, x=5/12 can be decom- 
posed as 1/2 © 1/2 © 1/3 © 1/2 © 1/3, where the © 
notation indicates a simultaneous addition of numerator 
and denominator in the fractions (e.g. n\/m\@ni/m,i = 
(rti +712) I (mi +^2), which physically corresponds to the 
concatenation of two subsets of motifs of concentrations 
xi=ni/mi and 2:2=712/7712) This leads to the so-called 
uniform sequence [22323] , which compactly indicates the 
composition and arrangement of the basic tiles in a pseu- 
douniform stripe configuration. 



reducible tiles of compositions 1/n and l/(n + 1) bracket- 
ing the overall B fraction xJ^. These tiles can be seen as 
encoding at a local level the drive for dilution, and their 
overall arrangement according to the rules results in a 
pseudouniform intergrowth. As discussed in Ref. [l7l the 
details of the repulsive model potential influence the de- 
tails of the stability of different arrangements, so a purely 
repulsive criterion does not univocally lead to a single so- 
lution in two dimensions. Hence our consistent use of the 
term pseudouniform to refer to the these orderings. Fine 
differences in stability are in practice irrelevant when one 
takes into account that the drive for dilution is a conve- 
nient but rough simplification of the state of affairs in 
real systems. This suggests that a fruitful approach to 
the problem of cataloging pseudouniform orderings in two 
dimensions is to abstract the basic ingredients, which are 
the irreducible, locally uniform, tiles, and their juxtaposi- 
tion into intergrowth arrangements following appropriate 
rules. 




FIG. 3: Stripe orderings on a square lattice of composition 
Ai- X B X for x = 5/12 (a) and x = 3/14 (b), obtained via 
simulated-annealing optimization using an isotropic Yukawa- 
type repulsive interaction among the B atoms.— (c) and (d): 
Fourier spectra (discs with radii proportional to the square 
modulus of the geometric structure factor) of the stripe ar- 
rangements shown in (a) and (b), showing the primary modu- 
lation wave vector. Only superstructure reflections are shown. 
The crosses indicate the reciprocal unit cell of the underlying 
square lattice. The superstructure reflections are indexed as 
satellite reflections according to the chosen primary modula- 
tion wave vector. 

Figure [3Jb) shows an analogous stripe ordering for 
x=3/14, where stripes with composition 1/5 and 1/4 
are arranged according to the uniform sequence [554]. 
The stripe pseudouniform sequence corresponding to any 
other composition can be derived using the Farey tree 
construction. In most cases the arrangement involves ir- 



Another key result of Rcf. [l3 is that the arrange- 
ments of the type shown in Figure [3] can be considered 
as compositionally modulated structures with a single 
composition-dependent primary modulation wave vector 
(a ID modulation)^ The Fourier spectra (square modu- 
lus of the geometric structure factor) of these two ex- 
amples are shown in Figure G3c) and EJd). One can 
see that the wave vectors (5/12)(l,l) and (3/14)(2,l), 
respectively, which arc approximately orthogonal to the 
stripe direction, determine the location of the most in- 
tense satellite reflections stemming from the ordering. 
Within the superspace formalism originally developed 
to describe incommensurately modulated structures ^r— 
these pseudouniform orderings can be described as long- 
period but commensurate structures with a simple step- 
like occupational modulation, with values either "atom 
A" or "atom 73" in the relevant proportion, and us- 
ing those vectors as primary modulation wave vectors. 
More specifically, a pseudouniform striped ordering is 
fully defined by an occupational modulation function 
f(xi) of period 1 with value "atom 73" in the inter- 
val — x/2 < X4 < x/2 and value "atom A" in the 
interval x/2 < X4 < 1 — x/2 and a modulation vector 
q that depends on x according to definite rules^ such 
that the kind of atom occupying a site m of the lattice 
is determined by the value of this function at a;4=q-m. 
The continuous variable x^ of the modulation function 
can be identified with the coordinate along the inter- 
nal space used in the superspace formalism. Thus, when 
described by means of occupational modulations, pseu- 
douniform orderings have their modulation functions re- 
duced to simple consecutive A- and 73-valued intervals 
(atomic domains) defined along the internal space in the 
superspace formalism^ 
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III. PSEUDOUNIFORM 2D COMPOSITIONAL 
ORDERING IN A PROTOTYPE Ai- X B X MODEL 

A. The snub-square ordering 

In some circumstances more complex, essentially two- 
dimensional, ordered configurations in the A\- X B X sys- 
tem can become competitive and prevail (in the sense 
of maximal dilution) over the ID striped arrangements 
studied in Ref. For the case of a square lattice and 
an ideal isotropic dilution drive this happens for composi- 
tions within the interval 1/5 < x <l/4. This composition 
interval is special because the basic tiles corresponding 
to the compositions x=l/A and x=l/5 can be juxtaposed 
with different orientations (the two sides of the two tiles 
have the same length), forming two-dimensional patterns 
(tilings). These 2D patterns were already considered in 
Ref. [22j in a simple context. Here we present a more com- 
plete picture of this family of orderings, their generating 
principles, and their description as modulated structures. 

Simulations for the repulsive lattice-gas model men- 
tioned in the previous section indeed confirm the appari- 
tion of intrinsically 2D ordered patterns in this composi- 
tion range. Fig. 2] shows the ground states obtained for a 
composition x~2/9 for two different sizes of the simula- 
tion supercell (subject to periodic boundary conditions). 
The pseudouniform [45] stripe arrangement, in agree- 
ment with the rules explained in the previous section, 
was obtained for a minimal 9x9 supercell (Figure |4jb)). 
But for a 12 x 12 supercell, a completely different ordering 
pattern is stabilized (Figure Eta)). It has the form of a 
so-called snub-square tiling^ (with a slight modification: 
the rhombic 1/4 tile can be considered formed by two 
triangles that are isosceles instead of being equilateral 
as in the canonical snub-square tessellation). This or- 
dering pattern (henceforth SSQ ordering) prevails in the 
simulations over the stripe arrangement if the boundary 
conditions are compatible with both types of orderings. 

A comparison of the distribution of B-B distances 
points to the main reason for the prevalence of the SSQ 
ordering over the striped one. As pointed out by Wat- 
son^ both have the same distribution of B-B inter- 
atomic distances up to \/l0 (in cell parameter units of the 
underlying square lattice), but the SSQ ordering avoids 
the next distance vT3 associated with neighbours sepa- 
rated by a vector (3,2) or (2,3), which is present in the 
striped arrangement. 

Figure @Jc) and E{d) presents sketches of the Fourier 
spectra of the SSQ ordering and of the competing striped 
arrangement. They clearly show the primary modula- 
tion wave vector(s) that can be identified in each case. 
The superstructure reflections of the striped arrangement 
can be indexed as linearly arranged satellites, with the 
strongest one being given by the modulation wave vector 
q= (2/9_U2,l), in accordance with the rules established 
in Ref. [l7|- The diffraction pattern of the SSQ ordering 
requires instead two primary wave vectors for a simple 
indexation of the satellite reflections, and the strongest 
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FIG. 4: Uniformity-driven orderings in a square lattice 
of composition Ai- X B X for x=2/9, obtained with simula- 
tions based on an isotropic Yukawa-type repulsive interaction 
among the minority B atoms, for 12x12 (a) and 9 x 9(b) sim- 
ulation boxes, and their corresponding Fourier spectra ((c) 
and (d)). The primary modulation wave vectors, one for the 
stripe arrangement (b) and two for the 2D tiling (a), are in- 
dicated in the corresponding Fourier diagram. The crosses in 
(c) and (d) indicate the reciprocal unit cell of the underlying 
square lattice. The superstructure reflections are indexed as 
satellites using the primary wave vector (s). 



satellites indicate the most obvious choice: qi = (l/2,l/6) 
and qa=(-l/6,l/2). 

The relevance of the SSQ-type arrangements as 2D or- 
derings that could maximize somehow the uniformity of 
the distribution of a minority component for some spe- 
cific compositions was pointed out in Ref. [22|. We stress 
here that these configurations arc indeed observed as 
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ground states in simulations of lattice gas models with re- 
pulsive interactions that mimick an isotropic drive to uni- 
formity. In the following we will show their hierarchical 
structure when they are interpreted in terms of uniform 
sequences, and how they can produce nano-chessboard 
and nano-diamond patterns. In this framework, it will 
be shown that this type of arrangements, when described 
as compositional modulated structures, have some sim- 
ple basic common ingredients which can be generalized 
and applied to explain and rationalize orderings in real 
systems generated at least in part by a drive to maximal 
uniformity, including ordered patterns in the nano scale. 



B. Generalization of the snub-square ordering for 
any composition between 1/4 and 1/5. 

The SSQ ordering for a;=2/9 can be seen as a perfect 
array of intercrossing ribbons, each one formed by single 
tiles arranged in a l/4©l/5 sequence. (From now on, 
to emphasize the actual arrangement of motifs in a se- 
quence, we will use a slightly modified notation: {45} in- 
dicates in this case an alternating sequence of 1/4 and 1/5 
tiles, while the [] form will be reserved for pure numeri- 
cal Farey term sequences for a specific concentration. In 
the stripe-ordering cases both notations are equivalent.) 
This topology, which avoids the distance vT3 present in 
the stripe arrangement, can be generalized to any other 
composition in the interval l/5< x <1/4i22 For compo- 
sitions different from 2/9, it is not possible to avoid the 
\/T3 distance completely, but in general the B-B distance 
distribution will still be favorable when compared to the 
stripe ordering. Following the argument in the previous 
section, we will discuss these generalized SSQ orderings 
in terms of our basic mechanism of juxtaposing appro- 
priate tiles following certain rules, and we will not be 
concerned any more with detailed questions of maximal 
stability within a repulsive model. Figure [5ja) shows 
an example of a generalized SSQ ordering with x <2/9. 
It corresponds to the case x=9/41, and is characterized 
by the crossing of two kinds of ribbons, one correspond- 
ing to the tile sequence {554} (shown in the figure with 
shaded 1/5 tiles) and the other to the sequence {445}. 
In each ribbon, the minority tiles appear in a propor- 
tion 1/3, and the {445} ribbons are also a minority in a 
proportion 1/3. The minority ratio 1/3 (which we will 
henceforth denote by a), is thus a key parameter in the 
construction of the pattern. In fact, if we take abstract 
motifs a and b, both the tile sequence in each ribbon 
and the ribbon arrangement sequence are instances of 
the physical sequence {aab} (uniform sequence [3]), as 
befits a concentration a=l/3. By comparison, for the 
SSQ ordering for £=2/9 in Figure E{a), a=l/2, which 
corresponds to a simple alternation of 1/4 and 1/5 tiles 
on each (identical) ribbon and a simple alternation in the 
ribbon arrangement (uniform sequence [2] for a=l/2). 

The parameter a fully characterizes the generalized 
SSQ orderings by determining the hierarchical tiling ar- 



rangement, both at the level of each ribbon, and at the 
level of ribbon kinds. As it relates to a minority concen- 
tration, we can take without loss of generality a <l/2, 
and the general relation between a and x is given by>^ 

X = 9±(l-2a) 2 (1) 

where the plus sign in the denominator corresponds to 
l/5< x <2/9 (with 1/4 tiles underrepresented) , and the 
minus sign to 2/9< x <l/4 (with relative abundance of 
1/4 tiles). 

Fig. [5{b) shows the corresponding arrangement 
for a composition much closer to 1/5, a=l/9 
(x=81/389=0. 20826). This ordering optimizes the dis- 
tribution of B-B interatomic distances towards larger 
values by forming large squares of composition a;=l/5, 
separated by single ribbons with an inverted 1/9 propor- 
tion of 1/5 tiles, in an effective sequence {444444445}. 
Despite the large size of the square blocks, the system 
cannot be considered subject to phase segregation. On 
the contrary, the regions with different local composi- 
tion are coherently interleaved and ordered as a single 
phase. Paradoxically, it is the quest for uniformity that 
can drive the system into this type of chessboard order- 
ings. For a values of type n/m, the patterns become more 
complex, as the ribbons with a proportion n/m of basic 
tiles of type either 1/4 or 1/5 must be ordered according 
to a ID non trivial pscudouniform sequence of the type 
explained in the previous section. The presence of ID 
uniform sequences (derived from the Farey-tree construc- 
tion) along the individual ribbons is in fact the signature 
that the ordering is being driven by a uniformity quest. 
Figure |5] shows two examples for a=2/7 and 3/7. One 
can see in this figure that a more complex uniform se- 
quence is realized in the way the 1/4 and 1/5 basic tiles 
are arranged along each oblique ribbon, with either 1/4 
or 1/5 tiles as minoritary in a proportion a. For a=2/7 
(=1/3© 1/4), the relevant abstract uniform sequence is 
{aabaaab} ([34]), leading to interlocking ribbons of the 
form {4454445} and {5545554}, which themselves are in 
a proportion 2/7 with the same sequencing. 

For a rational value of x 7 the corresponding a accord- 
ing to Eq. (1) is in most cases irrational, and therefore 
the generalized SSQ ordering for that x is incommensu- 
rate with respect to the underlying lattice. The ribbons 
of basic tiles will follow an aperiodic uniform sequence 
corresponding to the irrational value of a, still accord- 
ing to the Farey tree construction^ Thus, the tendency 
to uniformity could conceivably stabilize 2D incommen- 
surate orderings of minority motifs even if their propor- 
tion is a simple rational value. For instance, the gen- 
eralized SSQ ordering for x=3/14, which competes with 
the stripe arrangement shown in Figure EJb), would be 
an incommensurate ordering with a having an irrational 
value close to 2/9. 

Within the interval 2/9 < x <l/4 the number of 1/4 
tiles is larger than that of 1/5 tiles, and the ordering 
patterns follow then the same rules but with the roles of 




FIG. 6: Generalized snub square orderings in a square lattice 
of composition A X -«.B X for (a) 2=49/225=0.217778 (a=2/7) 
and (b) 2=49/221=0. 221719 (a=3/7) 



the intercrossing feature that tends to eliminate as much 
as possible the y/l3 distance. 



FIG. 5: Generalized snub square orderings in a square lattice 
of composition Ai~ x B x for (a) 2=9/41=0.219512 (q=1/3) 
and (b) 2=81/389=0.20826 (a=l/9) and their corresponding 
Fourier spectra ((c) and (d)). The crosses in (c) and (d) indi- 
cate the reciprocal unit cell of the underlying square lattice. 
The superstructure reflections are indexed as satellites using 
the primary wave vectors that are indicated. 

the two types of basic tiles interchanged, such that the 
ribbons having a proportion a of 1/5 tiles become ma- 
joritary. As x approaches the composition limit £=1/4, 
blocks of density 1/4 increase in size, forming patchwork 
patterns of diamond shape, as shown in Fig. [7J 

We see that the a parameter acts as an auxiliary "mi- 
nority concentration" , controlling the arrangement of the 
1/4 and 1/5 basic tiles, which in themselves encapsulate a 
maximal local uniform ordering. The only other essential 
ingredient that makes the SSQs orderings more uniform 
than the stripe arrangements for the same composition is 



C. Description of uniformity-driven 2D orderings 
as modulated structures 

Although the generalized SSQ orderings discussed 
above present a large variety of arrangements, their 
diffraction patterns are very similar. They exhibit a 
strong hierarchy in the intensity of the superstructure re- 
flections, characteristic of modulated structures (see Fig- 
uresHlini and[7|). The strongest superstructure reflections 
can be taken as first-order satellites, and define the pri- 
mary modulation vectors, while the remaining ones can 
be indexed as higher-order satellites. In Section MI Al 
we showed the primary modulation wave vectors that 
can be associated with the x=2/9 SSQ ordering. These 
are particular values among those that can be obtained 
for generalized SSQ arrangements, which depend on the 
composition according to simple rules of geometrical ori- 
gin. In FigurcJSJ;), for example, the second-order satellite 
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FIG. 7: (a) Generalized snub square orderings in a square lat- 
tice of composition Ai- X B X for 1=32/135=0.237037 (a=l/8, 
x >2/9) and (b) sketch of its Fourier spectrum. Superstruc- 
ture reflections are indexed as satellites with primary modu- 
lation vectors q'i and q'2 (which in this case, as x >2/9, do 
not follow Eqs. [2]) 



(11) lies along the line joining two first-order satellites, 
and at distances from them which are fractions a and 
1 — a of the length of the line. This constraint, together 
with a simple symmetry argument (existence of a sym- 
metry plane for x >2/9, and tetragonality for x <2/9), 
is enough to fix the modulation vectors. To be short, 
we only describe those corresponding to the cases with 
l/5< x <2/9. In this interval, the two modulation wave 
vectors are given by the equations: 



qi = ( 
q2 = (- 



1 



x(l 



2 

1 3a; 

2 + T- 



-2a) 1 3a; 
2 '2 _ ~2 
1 x(l - 2a) 



(2) 



and are an explicitly tetragonal and right-handed set, 
best suited for the construction of the superspacc model 
to follow. (There is an alternate choice of modulation 



vectors ki=qi and k2 =(1, 0) — q2, with the extra prop- 
erty that |k 1 xk 2 |=a;. Also, these "natural" vectors are 
reciprocal to the real space vectors denning a "virtual" 
monatomic unit cell that best approximates the uniform 
motif distribution with density x (in correspondence with 
the concept described in Ref. H3).) 

The SSQ ordering can now be described as a mod- 
ulated structure with the modulation vectors of Eqs. [5] 
and with an occupational modulation defined by a func- 
tion /(a;4,a;5), of period 1 for both variables, such that 
the atom occupying a given lattice site m=(mi,m2) is 
determined by the value (either "atom A" or "atom i?" ) 
of f(x4,x 5 ) for a;4=qi-m and a;5=q2-m. Maps of the 
function /(a^xs) sampled at the discrete values of X4 
and that correspond to the SSQ orderings considered 
so far are shown in Figure [5] As it happens in the ID 
case, the domain with value "atom B" , corresponding 
to the minority motifs, forms a single compact/dense 
region with an area equal to their proportion x in the 
lattice. Thus the minority atoms B aggregate in the in- 
ternal space of the occupational modulation function in 
order to describe arrangements in which these atoms are 
maximally scattered in real space. (As remarked above, 
for irrational a one would obtain from Eqs. [2] incom- 
mensurate wave vectors, and the sampling of the A- and 
-B-valued domains would be continuous.) The B atomic 
domains have in all cases a form which avoids the oc- 
cupation by B atoms of two neighbouring lattice sites. 
This restriction is ensured by the geometrical condition 
shown in Figure [8] The B atomic domains, if trans- 
lated on the plane (3:4, 25) by either (qi-(l,0),q2-(l,0)) 
or (qi-(0,l),q2-(0,l)), juxtapose with the original one, 
with no superposition. Within the superspace formal- 
ism this is termed "closeness condition" , and is known to 
be satisfied by the atomic domains describing some qua- 
sicrystalsi 24 ' 25 The atomic domain borders can be taken 
parallel to the basis vectors of the point lattice in inter- 
nal space formed by the points (qi-m,q 2 m), where m 
represents the direct lattice points. This, together with 
the closeness condition and Eqs. [2J is sufficient to define 
the appropiate atomic domain B for any composition x, 
and to construct with it its corresponding SSQ ordering. 



IV. UNIFORMITY-DRIVEN ORDERINGS IN 
REAL SYSTEMS 



We have considered above an idealized A\_ X B X model 
for which the optimal ordering within a lattice depends 
only on the maximization of the dilution of minority mo- 
tifs. In real systems this factor can be one among many 
others. For instance, intrinsic anisotropies of the un- 
derlying crystalline lattice, energy constraints on possi- 
ble orderings (for instance, compatibility with rigid unit 
modes, chemically forbidden geometries, etc.) can play 
a fundamental role in favoring a particular ordering pat- 
tern. Nevertheless, if an effective dilution still retains a 
main role in these realistic scenarios, some of the fea- 
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FIG. 8: Maps of the occupational step-like modulation func- 
tions f(x4,,X5) that describe, taking as modulation wave vec- 
tors the vectors qi and q2 of Eqs. [2] the generalized snub 
square orderings for x=2/9 (a), 9/41 (b), 81/389 (c) and 
49/221 (d), which are shown in Figures 14151 and [6] respec- 
tively. The (x4,xs) coordinates of B-occupation sites sites 
are indicated with larger dots, and A-sites with smaller ones. 
The A (dark) and B-valued (clear) domains of the simplest 
function consistent with the occupied A and B sites are in- 
dicated in each case. A single unit cell of the internal space 
is shown. The closeness condition satisfied by the neighbor- 
ing B-valued domains in the superspace description is shown 
by means of continuous lines corresponding to their projected 
borders. 



tuxes of the purely uniformity-driven orderings are very 
likely to be maintained. One of the features which are 
present in the above discussion is the possibility of con- 
structing the orderings through the juxtaposition (inter- 
growth) of basic tiles. These encapsulate at the local 
level the drive for dilution, and their combination, as we 
have seen, can lead to a variety of patterns, depending 
on the geometrical flexibility involved. Fig. Ela) shows 
a sketch of the ordering pattern of oxygen vacancies ob- 
served in layers of the compounds Las-zSr^CusC^o-i, 
(La,Sr) 8 CugOi8 and La2Sr 6 Cu80i6^ To our knowledge, 
this is the first time that this ordering is interpreted as a 
tiling, similar to the SSQ, but with tiles of composition 
1/3 and 1/5, so that the global vacancy density is 1/4. In 
this case the ordering is not optimizing a purely repul- 
sive interaction (which would lead for this composition 
to a monatomic superlattice of vacancies), but the drive 
for dilution is still obviously present. Considered as an 
abstract geometrical and numerical exercise, the combi- 
nation rules for the 1/3 and 1/5 tiles lead to the same 
kind of regularities discussed above for the 1/5, 1/4 tiles, 
including the appearance of coherent generalized nano- 
patches (nano-rhombi in this case)i2I Crucially, all the 
arrangements produced through the 1/5-1/3 tile combi- 
nations are also describable as occupationally modulated 
structures, with compact atomic domains in the super- 
space construction, and with modulation vectors that can 



be read directly from the diffraction diagram. For the 
case of the vacancy distribution in Fig.j9ja), the primary 
modulation wave vectors can be chosen as qi=(l/2,l/4) 
and q2=( — 1/4, 1/2), as seen in the diffraction diagram 
of Fig. |9jb) , and the corresponding modulation function 
is shown in Figure EJc). The four independent vacan- 
cies within the 16 atomic sites form again a single simple 
atomic domain fulfilling the closeness condition described 
above. 

Once again we see that the ordered patterns are de- 
scribable as modulated phases with step-like modulation 
functions, with the £?-valued regions of the functions lim- 
ited to one or a few domains within the unit-cell of the 
periodic function, and with the fulfilment of some kind 
of closeness condition. Modulated structures with mod- 
ulation functions of this type arc very efficient for dis- 
tributing regularly minority motifs in real space. With 
all the evidence presented so far (here and in the previous 
work over several years )2i2§ri£ we can say that in a real 
system the knowledge of the relevant modulation vectors 
(extracted from its diffraction diagram) is sufficient to de- 
sign, through the rules explained above, a composition- 
ally modulated model which corresponds to an atomic 
ordering with a high degree of uniformity in the distri- 
bution of some specific motif. Thus we can postulate 
an atomic distribution model with a specific composi- 
tion and a maximal dilution consistent with the observed 
modulation vectors. If a drive for dilution is at work, this 
apriori model should approximate the experimental ar- 
rangement. This clearly occurs in the example shown in 
Fig. HI and also seems to happen in the nano-chessboard 
arrangements observed in the compounds mentioned in 
the introduction. As shown below, the basic features 
of the latter can be reproduced in a simplified binary 
system by postulating a compositionally modulated ar- 
rangement with step-like atomic domains that fulfill the 
closeness condition, and are consistent with the observed 
modulation vectors. This would indicate that the ob- 
served arrangements represent some optimization of the 
dilution of the atoms within the underlying perovskite 
lattice, but subject to some anisotropic contraints that 
force the possible modulation vectors. 

Let us consider for instance the compound 
Lii/2-3i/Nd 1 /2+yTi03 mentioned in the introduc- 
tion. Figure [TU] shows a [001] zone electron diffraction 
pattern (EDP) of this material with y=0.067 (we change 
the letter for the composition variable to avoid confusion 
with that employed in our binary toy model). The 
primary modulation wave vectors can be identified 
directly from the EDP diagram as qi=(l/2— e,— 1/2) and 
q2=(l/2,l/2— e) with e about 1/30. One cannot pretend 
to derive from this information a quantitative structural 
modelling of this compound, but we can determine for 
our simple 2D A\- X B X binary system a pseudouniform 
ordering of the B atoms consistent with these modulation 
vectors. Using the rules explained above, we have to 
consider an occupational modulation function f(x^,x§) 
fulfilling the closeness condition for these modulation 
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FIG. 9: (a) Ordering pattern of the oxygen vacan- 
cies reported in the layers BO2-S of Lag-xSr^CusCbo-a , 
(La,Sr)gCu80ig, and La2Sr6CugOi6 (Ref. [2r3 ). interpreted 
as a 1/301/5 tiling in a virtual square lattice in which va- 
cancies are present in a proportion x=l/A. (b) Diffraction 
diagram showing the modulation vectors qi = (1/2, 1/4) and 
q2=(-l/4,l/2) that index the satellites, (c) Map of the mod- 
ulation function that can be associated with the ordering pat- 
tern using qi and q2 (symbols and extra closeness-condition 
lines as in Figure [HJ. 



of this A\_ X B X nano-chcssboard arrangement is fully 
determined by the modulation vectors, and is given by 
the relative area of the B atomic domain in Fig. fTTT b). 
which is x=l/2-e + e 2 , i.e. x=421/900=0.467778. 
Obviously this composition is not comparable with that 
of Lii/2-3yNdi/2+j,Ti03. The arrangement of Fig. QT] 
is far from the real system, not only due to the reduc- 
tion to a 2D Ai_ x B x arrangement, but also because 
the vacancies accompanying the Nd/Li substitution 
and the expected very large positional relaxations are 
ignored. But nevertheless, this simple occupational 
modulation with an intrinsic tendency to uniformity 
and consistent with its modulation vectors in a binary 
system, is sufficient to reproduce basic features of the 
real system. Notice for instance the shift of the A and 
B occupation sites in contiguous patches of composition 
^4i/2-Si/2 of the chessboard, as proposed in the model 
for Li 1/2 _ 3 . y Nd 1/2+ . y Ti0 3 of Ref. Q. Fig. QTIc) shows 
the geometric diffraction pattern of this 2D A\_ X B X 
nano-chessboard arrangement. It is remarkable that 
despite the strong anharmonicity of the occupational 
modulation only odd-order satellites are observable. 
This is consistent with the Fourier decomposition of 
the underlying two dimensional step-like occupational 
modulation. In the real system even-order satellites 
close to the main reflections are also significant, which 
is probably due to the strong displacive modulations in 
the perovskite framework (mainly tiltings), induced by 
the cation ordering. This can be a plausible explanation, 
since by just introducing a small displacive sinusoidal 
modulation with wave vectors qi+q2 and qi-q 2 of the 
A and B positions in the configuration of Fig. [TTT a). 
we can produce a diffraction diagram similar to the 
experimental one of Figure 1101 

Nano-chessboards of the type reported in Refs. [1|- 
may therefore originate in a drive for maximal unifor- 
mity in the distribution of a minority motif. If this were 
indeed the mechanism, the width of the interface regions 
between the chessboard patches, although clearly wider 
than in our model due to structural displacive relax- 
ations, should be independent of the size of the patches, 
i.e. of the system composition. A careful check of the ful- 
fillment or not of this property by the nano-chessboards 
observed in real systems would be a key feature for the 
validation of a uniformity-drive mechanism. 



V. CONCLUSIONS 



vectors. This simple step-like occupational modulation 
function (see Fig. lllf b)). yields the atom ordering shown 
in Fig. [TTta). which vividly recalls the nano-chessboard 
arrangements observed in this compound. The square 
patches have composition AB, while the global excess of 
A atoms are localized at the interfaces, which are made 
only of A atoms and are limited to two unit cells of the 
underlying lattice. By construction, the composition 



We have shown that a drive towards dilution of minor- 
ity structural motifs constrained to an underlying lattice 
is sufficient to produce two-dimensional orderings with 
local composition variations in the nanoscale. The di- 
lution can be achieved in practice by the combination 
or intcrgrowth of basic tiles that encapsulate a local 
repulsion of minority motifs. This combination princi- 
ple is able to generate orderings that are hierarchically 
structured, adopting for certain compositions the form of 
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FIG. 10: [001] zone electron diffraction pattern of the (3+2)- 
D incommensurately modulated, nano chessboard phase of 
Lio.3oNdo.567Ti03, indexed with respect to the basis set M* 
= { a *= a p, b * = b*, c*= l/2c*, qi = l/2a* + (l/2-e)bp , 
q2 = (-l/2+e)ap + l/2bp}, where the subscript p stands for 
parent perovskite, and e ~ 1/30. 
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nano-chessboards or diamond-like patches. In all cases, 
regions of two different effective compositions are inter- 
leaved coherently, with long-range order. 

The mechanism leading to composition patterning 
from a dilution drive has already been documented in 
other systems: layer arrangement in compositionally flex- 
ible layered compounds^i^Sr— and two-dimensional sys- 
tems exhibiting effectively ID stripe patterns^ Here we 
find what appears to be a segregation of phases instead 
of a pscudouniform distribution, but these arrangements 
arc in fact ncar-optimally uniform by any reasonable cri- 
terion, and their underlying structure (now essentially 
two-dimensional instead of effectively ID) can be de- 
scribed by the same recipe at work in the pseudouniform 
systems previously studied: modulation vectors directly 
and simply determined from the intensity distribution in 
the diffraction diagram, and step-like occupational mod- 
ulations satisfying a closeness condition for these modu- 
lation vectors. 

The description in terms of occupational modulations, 
best represented using the tools of the superspacc formal- 
ism, is thus seen as the unifying structural principle in a 
wide variety of systems which exhibit in some degree a 
dilution drive. The key idea is that minority motifs max- 
imally scattered in real space are represented by compact 
regions in the internal coordinates of the superspace con- 
struction. 

We have shown that this powerful principle seems 
to be at work for two kinds of relevant experimen- 
tal systems. In the case of the distribution of 
vacancies in Lag-xSr^CugC^o-s, (La,Sr)gCu80i8 and 
La2Sr@Cu80i6r^ the resulting pattern can be seen also 
as a direct example of the tiling rules in real space. 




FIG. 11: (a) Chessboard nanopattern in a square lattice of 
composition A\- X B X , obtained as a compositional modulated 
structure with modulation wave vectors qi=(l/2,l/2— e) 
and q2 = (-l/2+e,l/2) with e=l/30 (i.e. those present in 
Lii/2-3»Ndi/2+ y Ti03), and a step-like A/B occupational 
modulation function fully determined by the closeness con- 
dition of the B domains for these modulation vectors. This 
modulation function is shown in (b), where A and B-valued 
domains are depicted as dark and clear regions, respectively, 
within a 2D period of the function (i.e. a unit cell of the inter- 
nal space in the superspace construction) . The closeness con- 
dition satisfied by the neighboring B-valued domains in the 
superspace description is shown by means of continuous lines 
corresponding to their projected borders. The corresponding 
geometric diffraction pattern is shown in (c). 



The nano-chessboard arrays, and similar patterns ex- 
hibiting composition patterning, that have been re- 
cently observed in composition flexible systems such as 
Lii/2-3a;Nd 1 /2+z:Ti03 are strongly reminiscent of the 
kinds of hierarchical real-space orderings that appear in 
uniformity-driven simple models. They reproduce basic 
features of the observed nano-chessboard arrangements, 
and their Fourier spectra exhibit the kind of strong hi- 
erarchical satellite structure present in the experimental 
diffraction diagrams. We can thus suggest that the ob- 
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served compositional changes at the nanoscale in these 
compounds have their origin, at least partially, in a drive 
for uniformity, maintaining the coherency of a single 
phase, and are not, as previously proposed, due to any 
kind of phase separation. 
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